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The nature of finite temperature transitions in lattice QCD with Wilson 
quarks is studied near the chiral limit for the cases of 2, 3, and 6 flavors of 
degenerate quarks {Np = 2, 3, and 6) and also for the case of massless up and 
down quarks and a light strange quark [Np = 2 + 1). Our simulations mainly 
performed on lattices with the temporal direction extension A^t = 4 indicate 
that the finite temperature transition in the chiral limit (chiral transition) 
is continuous for Np = 2, while it is of first order for Np = 3 and 6. We 
find that the transition is of first order for the case of massless up and down 
quarks and the physical strange quark where we obtain a value of m^/nip 
consistent with the physical value. We also discuss the phase structure at 
zero temperature as well as that at finite temperatures. 



1 Introduction 



One of major goals of numerical studies in lattice QCD is to determine the 
nature of the transition from the high temperature quark-gluon-plasma phase 
to the low temperature hadron phase, which is supposed to occur at the early 
stage of the Universe and possibly at heavy ion collisions. It is, in particular, 
crucial to know whether the transition is a first order phase transition or a 
smooth transition (second order phase transition or crossover) to understand 
the evolution of the Universe. 

Determination of the order of the transition for the case of degenerate 
Np flavors, is an important step toward the understanding of the nature of 
the QCD transition in the real world. We can compare the numerical results 
for various number of flavors with theoretical predictions based on the study 
of the effective a model |l], 0. In order to investigate what really happens 
in the nature, we have to ultimately study the effect of the strange quark 
together with those of almost massless up and down quarks, because the 
critical temperature is of the same order of magnitude as the strange quark 
mass. 

In this article we investigate finite temperature transitions in lattice QCD 
using the Wilson formalism for quarks for various numbers of flavors {Np = 2, 
3, and 6) near the chiral limit and also for the case of massless up and down 
quarks and a light strange quark {Np = 2 + 1). Most simulations of finite 
temperature QCD were performed with staggered quarks. However, because 
the Wilson formalism of fermions on the lattice is the only known formalism 
which possesses a local action for any number of flavors, it is important to 
investigate the finite temperature transition with Wilson quarks and compare 
the results with those for staggered quarks. 

In Sec. H, we define our action and coupling parameters. Because chiral 
symmetry is explicitly broken on the lattice in the Wilson formalism, we 
first define the chiral limit for Wilson quarks and give a brief survey of the 
phase structure in Sec. ^ Our simulation parameters are summarized in 
Sec. |[ Numerical results for the chiral limit are summarized in Sec. |^. We 
then discuss, in Sec. ^ problems and caveats which appear in a study of 
the finite temperature transition with Wilson quarks when performed on 
lattices available with the present power of computers. Sec. |^ deals with 
the transition in the chiral limit (chiral transition) in the degenerate cases 
of Np = 2, 3 and 6. In Sec. ||, we study the influence of the strange quark 
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on the QCD transition both in the degenerate Np = 3 case and in a more 
realistic case of the massless up and down quarks with a massive strange 
quark, Np = 2 + 1. We finally conclude in Sec. ^. Preliminary reports are 
given in P, I |. 

2 Action and coupling parameters 

We use the standard one-plaquette gauge action 

S, = ^j:^eTTiUp) (1) 
y p 

and the Wilson quark action ^ 

Np 

-5, = - E E M^) D{Kf, n, m) ^^(m), (2) 

f=l n,m 

D{K,n,m) = Sn,m-K^{{l - 7m) f^n.^f^n+^.m + (1 + 7m) ^m,/.'^m+M,n}> (3) 

where g is the bare coupling constant and K is the hopping parameter. In 
the case of degenerate Np flavors, lattice QCD contains two parameters: the 
gauge coupling constant /? = and the hopping parameter K. In the 
non-degenerate case, the number of the hopping parameters is Np. 

We denote the linear extension of a lattice in the temporal direction by 
Nt and the lattice spacing by a. 

3 Brief survey of phase structure 

In the Wilson formalism of fermions on the lattice, chiral symmetry is ex- 
plicitly broken by the Wilson term even for vanishing bare quark mass [^]. 
The lack of chiral symmetry of chiral symmetry causes much conceptual and 
technical difficulties in numerical simulations and physics interpretation of 
data. Therefore before going into discussion of details of data and analyses, 
we give a brief survey of the phase structure at zero temperature as well 
as that at finite temperatures |0, |], including the results presented in this 
article. 
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3.1 Quark mass and PC AC relation 

We first define the quark mass through an axial- vector Ward identity 0, p!0| ]. 

2m, (0|P|7r(p = 0)) = -m^ ( | 1 7r(p = 0)) (4) 

where P is the pseudoscalar density and ^4 the fourth component of the 
local axial vector current. [Note that we have absorbed a multiplicative 
normalization factor into the definition of the quark mass m,, because this 
convention is sufficient for our later study. We also note that there is an 
alternative definition of the quark mass replacing with e.g. (exp(— m^a) — 
l)/a, which gives the quark mass identical with the above within order of a.] 
With this definition of quark mass, the PCAC relation. 
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(5) 



which is expected to be satisfied near the continuum limit, was numerically 
first verified within numerical uncertainties for the quenched QCD at zero 
10], |lT| and subsequently for various cases including QCD 
i. 111, [T|, H, |l6l. It should be noted that the PCAC 
relation is satisfied not only in the continuum limit, (3 = 00, but also even in 
the strong coupling limit, /? = 0: The result of the strong coupling expansion 



temperature in 
with Np = 2 in 



without quark loops 



cosh(m7ra) 



2mqa 



m-n-a ■ 



[1 - 16K^){1 ~ 4K^) 
AK^ (2 - 12^2) 
AK'^ sinh(m7ra) 
— 4i^2 cosh(m,ro) ' 



(6) 



gives the relation oc uiq at small m,. Our numerical data for Np = 2 at 
/3 = agrees well with these formulae within errors as shown in Fig. 
We note that, if 

(0|A4|7r(p = 0)) ocm^ (7) 

is satisfied for small as is the case both for /5 = and /? = 00, then the 
definition implies that the PCAC relation (|) is exact. It should be also 
noted that Eq.(|^ holds when Euclidean invariance is recovered [|10 . 



^In Ref . [p2[ , agreement between Eq.(^ and numerical data in the confining phase is 
shown also for the case Np = 18. The rho meson mass, the nucleon mass, and the delta 
mass also agree with corresponding strong coupling mass formulae. 
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Eq.(^ implies that when = 0, either = or ( | | 7t{p = 0) ) = 0. 
This further imphes, when we define the pion decay constant by 

(0|A4|7r(p = 0)) =m^/^, (8) 

that when = 0, either = or = is satisfied. Note that /,r = is 
the relation which should be satisfied when chiral symmetry is restored, and 
that = is the relation when chiral symmetry is spontaneously broken, 
both in the chiral limit. It might be emphasized that although the action 
does not possess chiral symmetry, either relation of m,r = or /^r = holds 
in the massless quark limit when the quark mass is defined by Eq.(^). In 
particular, in the confining phase, = when niq = and vice versa. 

3.2 Definition of chiral limit and phase structure at 
zero temperature 

We identify the chiral limit as the limit where the quark mass vanishes at 
zero temperature. This defines a chiral limit line Kc in the K) plane, 
which is a curve from /r~l/4at/5 = 0to-ft' = l/8at/5 = oo. See Fig. |[ 
In the following we also discuss alternative identifications of the chiral limit. 
When clear specification is required, we denote this Kc as Kc{mq). 

Let us denote a line where the pion mass vanishes at zero temperature 
by Kf.{m'^). This line is the critical line of the theory because the partition 
function has singularities there. As discussed in the previous subsection, 
we expect that Kc{mg) and Kc{m'^) are identical for small Np. It should 
be, however, noted that the Kc{mq) line is conceptually different from the 
Kc{m\) line: If quarks are not confined and chiral symmetry is not sponta- 
neously broken, there is no Kc{m1) line. In fact, for the case of Np > 7, the 
Kc{mq) line belongs to the deconfining phase and remains nonzero there 
— i.e. there is no Kc{m1) line around the Kc{mq) line, at least for small (3 

As a statistical system on the lattice, QCD with Wilson quarks is well- 
defined also in the region above the Kc line. Some time ago, S. Aoki 
proposed and numerically verified that the critical line Kc{m1) (for small Np) 
can be interpreted as a second order phase transition line between the parity 
conserving phase and a parity violating phase. This interpretation is useful 
in understanding the existence of singularities of the partition function. Once 
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its existence is established, various properties of hadrons can be investigated 
in the parity conserving phase. In particular, even with the Wilson term, 
various amplitudes near the chiral limit do satisfy Ward-Takahashi identi- 
ties derived from chiral symmetry to the corrections of 0{a) Therefore, 
although the action does not have chiral symmetry, the concept of the spon- 
taneous breakdown of chiral symmetry is phenomenologically very useful. 
Because our main interest is to study the physical properties of hadrons in 
the continuum limit, it is important to study these axial Ward-Takahashi 
identities and estimate the magnitude of the 0(a) corrections from the Wil- 
son term in the physical quantities. 

We have defined the Kc line by the vanishing point of at zero tem- 
perature, because this hne corresponds to massless QCD. In this connection, 
however, it should be noted that there necessarily are ambiguities of 0(a) off 
the continuum limit for lines in the (/?, K) plane which give the same theory 
in the continuum limit. This is true also for massless QCD: Instead of the 
condition m^r = 0, we may fix other quantities such as mp/rriN, which will 
lead to a line different from the Kc line. Of course, the continuum limit is not 
affected by these 0(a) ambiguities. We, however, would like to stress that 
the definition we have taken for the Kc is conceptually natural and useful for 



the reasons given in Sec. 3.1 



3.3 Phase structure at finite temperatures 

The temperature on a lattice with the linear extension in the temporal di- 
rection Nt is given by T = l/Nfa. On a lattice with a fixed Nf, finite 
temperature transition or crossover from the low temperature regime to the 
high temperature regime occurs at some hopping parameter when (3 is fixed. 
This defines a curve Kt in the (/3, K) plane. In this paper, for simplicity, 
we use the term "transition" for both genuine phase transitions and sharp 
crossovers, unless explicitly specified. At finite temperatures we denote the 
screening pion mass by m,r and sometimes we call it simply the pion mass, 
and similarly for other hadron screening masses. Quark mass at finite tem- 
peratures is defined through Eq.(^ with m.^ the screening pion mass, and 
similarly for through Eq.(^). Note that, with these definitions of m^r and 

^ In the particular form of Eq.(Q), we have absorbed these 0(a) corrections in the 
definition of m^, or, equivalcntly, in the value of Kc- 
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/tt, the discussions given in Sec. hold also at finite temperatures. 

One of fundamental problems is whether the finite temperature transition 
line Kt does cross the chiral limit line K^, where we define the K^. line by 
the vanishing point of iriq at zero temperature (cf. Sec. p.2|) . If the Kt line 
does not cross the Kc line, it means that there is no chiral limit in the low 
temperature confining phase. Therefore it is natural to expect that it does 



cross. However, as first noted by Fukugita et al. |]18| , it is not easy to confirm 
this: The Kt line creeps deep into the strong coupling region. In this paper 
we show that the Kt line indeed crosses the chiral line ii'cat/3~3.9 — 4.0 
at A^t = 4 and /5 ~ 4.0 — 4.2 at A^t = 6 for the case of Np = 2. (For previous 
reports see Refs. P, ||.) 

Because the Kc line describes the massless QCD, we identify the crossing 
point of the Kc and Kt lines as the point of the finite temperature transition 
of the massless QCD, i.e. the chiral transition point. (We will discuss later 
0{a) ambiguities in the definition of the chiral limit at finite temperatures 
which come from the lack of chiral symmetry.) 

Numerical studies show that, in the confining phase, the pion mass van- 
ishes, for a fixed (3, at the hopping parameter which approximately equals 
the chiral limit Kc. On the other hand, in the deconfining phase, the pion 
mass is of order of twice the lowest Matsubara frequency 271 /Nt in the chiral 
limit. Therefore, in the deconfining phase, the system is not singular even 
on the Kc line. 

Recently, Aoki et al. ||19|| investigated a critical line where the screening 
pion mass vanishes at finite temperatures, which we denote by Kc{m\] T ^ 
0). Based on analytic studies of the 2d Gross-Neveu model and numerical 
results in lattice QCD with Np = 2, they showed that the Kcirn^; T 7^ 0) 
line starting from i^'~l/4at/5 = sharply turns back upwards (to larger 
K region) at finite j3. The lower part of the Kc{m'^]T ^ 0) line is almost 
identical with the Kc{m\) line up to the sharp turning point, while the 
analytic results of the 2d Gross-Neveu model suggest that they slightly differ 
from each other, probably with 0(a). See Fig. 0. 

The non-existence of the Kc{m'^', T 7^ 0) line in the large j3 region is con- 
sistent with the previous results that m^r does not vanish in the deconfining 
phase along the chiral line Kc. The slight shift of the Kc{m1; T 7^ 0) line 
from the Kc{m'^) line in the confining phase was observed also in our pre- 
vious study [D, 1^ (see also Sec. |). This slight shift of the Kc{rn\]T 7^ 0) 
line means that is not rigorously zero on the Kc{m\) line in the confining 
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phase at finite temperatures. This small pion mass on the Kc line in the 
confining phase is caused by the chiral symmetry violation due to the Wilson 
term and should be of 0(a). 

Similarly to the Kc{m'^;T ^ 0) line, we define the line Kc{mg;T ^ 0) 
where the quark mass vanishes at finite temperatures. When we follow the 
line Kc{mq\T ^ 0) from /5 = 0, it is first identical with the Kc{ml.]T ^ 
0) line. The line Kc{mq,T ^ 0) passes through the turning point of the 
Kclm^; T 7^ 0) line and runs into the larger jS region, where /^r starts to vanish 
instead of on the Kdrnq, T 7^ 0) line. See Fig. ||. This suggests that the 
turning point which is the boundary between = and = is the finite 
temperature transition point. This further implies that the transition line 
Kt touches the turning point of the Kdrn^] T 7^ 0) line and moves upwards 
in the (/?, K) plane. This observation is not in accord with the argument by 
Aoki et al. |]T^ that there is a small gap between the K(.{m'^\ T 7^ 0) and Kt 
lines. 

We have identified the crossing point of the and Kt lines as the chiral 
transition point. In connection with the 0(a) ambiguities of the line for 
massless QCD in the coupling parameter space mentioned in Sec. p.2| , there 
are 0(a) ambiguities also in the definition of the chiral transition. Therefore, 
one may alternatively identify the sharp turning point of the Kc{'m^] T 7^ 0) 
line as the chiral transition point. The property of the chiral transition in 
the continuum limit is, of course, not affected by these 0(a) ambiguities. 



3.4 Characteristics for Wilson quarks 

Let us summarize several characteristic properties of the phase diagram of 
QCD which are originated from the explicit chiral symmetry violation of the 
Wilson term. They are in sharp contrast with those of staggered quarks 
where at least a part of chiral symmetry is preserved. 

(i) In the coupling parameter space, the location of the point where 777,7^ = 
in the confining phase is not protected by chiral symmetry off the continuum 
limit. Therefore, the chiral limit Kc-, defined by iriq = or m^r = at zero 
temperature, is different from the bare massless limit = 1/8 except at 
(3 = 00. 

(ii) As a statistical system on the lattice, QCD with Wilson quarks is well- 
defined also in the region above the K^. line. At zero temperature, the K^ line 
is a second order transition line between the conventional parity conserving 
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phase at K < Kc and a parity violating phase at K > Kc ^ 

(iii) At finite temperatures, the critical line Kcirn^; T 7^ 0) where the 
screening pion mass vanishes is not a line from ~ 1/4 at /? = to an end 
at some finite (3, but it sharply turns back toward larger K region at the 
finite /3 [||]. 

(iv) Although the major part of the effects from the Wilson term can 
be absorbed by the shift of from K = 1/8, there still exist additional 
small 0{a) effects which are related to the chiral symmetry violation. In 
particular, the location of the point where m^r = in the confining phase 
slightly depends on Nt |]T9[. The continuum limit is not affected by these 
0(a) effects. 



4 Simulation Parameters 

In this article we mainly perform simulations on lattices with the temporal 
direction extension Nt = 4. The spatial sizes are 8^ x 10 and 12^. To 
study the Nt dependence for the Np = 2 case, we also make simulations 
on Nt = 6 and 8 lattices. Simulations on an A^^ = 8 lattice are performed 
also for the case of Np = 2 + 1. When the hadron spectrum is calculated, 
the lattice is duplicated in a direction of lattice size 10 or 12. We use an 
anti-periodic boundary condition for quarks in the t direction and periodic 
boundary conditions otherwise. 

We generate gauge configurations for Np = 2 by the Hybrid Monte Cairo 
(HMC) algorithm with a molecular dynamics time step At chosen in such 
a way that the acceptance rate is about 80 — 90%. For Np > 3 and Np = 



2 + 1 we use the hybrid R algorithm with At = 0.01, unless otherwise 



stated. We fix the time length of each molecular dynamics evolution to r = 1. 
The R algorithm introduces errors of O(Ar^), while the HMC algorithm is 
exact. As reported recently also for staggered quarks [^, we note that step 



size errors with the R algorithm are large in the confining phase near the 
chiral limit. In the immediate vicinity of the chiral transition, we observe 
step size errors also in the deconfining phase where a large At can even push 
the phase into the confining phase, as reported previously with staggered 
quarks [^. In these cases, we apply a sufficiently small At so that the 
results for physical quantities become stable for a change of At. 

The inversion of the quark matrix is done by the minimal conjugate resid- 
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ual (CR) method with the ILU preconditioning or the conjugate gradient 
(CG) method without preconditioning. We find that the CR method is ef- 
ficient in the confining phase when it is not very close to the chiral hmit 
and also in the deconfining phase at large j3 and small K. In other cases we 
use the CG method. The convergence condition for the norm of the residual 



r is Y ||r|p/(12V^) < 4.5 x 10 (1.0 x 10 ^) for configuration generations 
(hadron measurements), where V is the lattice volume. We also check that 
the relative changes of the quark propagator at several test points on the x 
and t axes are smaller than 10"'^ for the last iteration of the matrix inversion 
steps: \{Gn — Gn-i)/Gn\ < 10"'^ whcrc n denotes the last iteration. In the 
HMC calculations, we check that the difference of the action after molecular 
dynamic evolutions is sufficiently small with this convergence condition. 

The statistics is in general totally r = several hundreds. The initial con- 
figuration is taken from a thermalized one at similar simulation parameters 
when such a configuration is available. In most cases, the plaquette and the 
Polyakov loop are measured every simulation time unit and hadron spectrum 
is calculated every 5t = 10 (or less depending on the total statistics). When 



the value of /5 is small the fiuctuation of physical quantities are small [|I2 
and therefore we think the lattice sizes and the statistics are sufficient for 
our purpose to determine the global phase structure of QCD at finite tem- 
perature. Errors are estimated by the single-elimination jackknife method. 
Simulation parameters are summarized in Tables |1| — ^ 



5 Numerical results for Kc 

As discussed in Sec. |3.2| , the chiral limit Kc is defined by the vanishing point of 
niq at zero temperature. One straightforward way to determine numerically 
the chiral limit at a fixed value of (3 is to calculate the the quark mass through 
Eq.(0) at several hopping parameters and extrapolate them to its vanishing 
point in terms of a linear function ofl/K. We denote the thus determined 
by Kcijnq). Because we expect the PCAC relation (|^) to hold also at finite 
/3, we may alternatively calculate K^. by the vanishing point of using a 
linear extrapolation of in 1/K. We denote this K^. by Kc{m'^). 

On finite temperature lattices, it was previously shown that the value 
of the quark mass at given {j3, K) does not depend on whether the system 
is in the deconfining or confining phase at (3 = 5.85 in the quenched QCD 
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1^ and at 13 = 5.5 for the Np = 2 case |T3[. This enables us to determine 
the chiral hmit, for these values of f3, alternatively by the vanishing point 
of rriq at finite temperatures. Strictly speaking there are systematic errors 
which come from finite Nf, as mentioned earlier. On the other hand, in the 
deconfining phase, one is able to perform simulations around the line as 
discussed later, i.e. we can determine without an extrapolation which 
usually leads to a considerable amount of systematic errors. Therefore, the 
determination of Kc from in the deconfining phase is useful in particular 
at large p. 

At small f3 region {(3 < 5.3) where we mainly perform simulations in this 
work, nig in the deconfining phase does not agree with that in the confining 
phase. Therefore, the proportionality between ruq in the deconfining phase 
and ml in the confining phase is lost, contrary to the case /? > 5.5 discussed 
above. This behavior is seen in Figs. ^ and ^ where physical quantities for 
Np = 2 at /3 = 5.0 and 4.5, respectively, are shown. As we discuss in Sec. ^ 
we interpret this unexpected phenomenon at /3 < 5.3 in the deconfining phase 
as a lattice artifact. 

In the confining phase, on the other hand, the proportionality between rriq 
and ml is well satisfied for all values of P 0, |12|, 0, |T6[ . We also find 



that mq and m^^ are almost independent of Nt in the confining phase. See 
Fig. ^ for Np = 2 at /3 = 4.5. Therefore we can calculate Kc approximately 
also by the vanishing point of mq, Kc{mq), or that of m^, Kc{m1), in the 
confining phase at T > 0. 

The numerical results for for Np = 2 obtained by various groups 
[|, 1^, 1^, |2^ are plotted in Fig. ^ together with finite temperature 

transition lines discussed in the following sections. The values of Kc show 
a slight dependence (at most of the order of 0.01) on the choice of Kc{mq) 
or Kc{ml), which can be probably attributed to the systematic errors in the 
extrapolation of and mq in l/i^',0 because, as discussed above, we expect 
that Kc{mq) and Kc{m1) are identical. The values of Kdml) for A''^;' = 2 for 



various /9's are listed in Table |I0| We estimate the systematic errors due to 
the extrapolation are of the same order as the differences between Kc{mq) 
and Kciml). 



The range of the quark mass value we use in this article for the extrapolation to 
determine the Kc is mainly about 0.2 — 0.5 in lattice units in the confining phase. As 
seen from Fig. ^, and vriq sometimes show slightly convex curves in 1/ K . In such cases, 
a choice of the fit range at smaller rUq will lead to slightly smaller values for K^. 
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The Nt dependence of Kc at P = 4.5 are listed in Table [Tl|. The Np 
dependence is also given. We find that the differences due to Np and Nt 
are of the same order of magnitude as the difference between Kc{mg) and 

To summarize this section, we note that although the chiral limit is de- 
fined by the vanishing point of rriq at zero temperature, there are several 
practically useful ways to determine Kc- Kdniq) and Kdm^) at T = and 
in the confining phase, and K(.{mq) in the deconfining phase. They all gives 
the same results within present numerical errors. 



6 Finite Temperature Transition and Prob- 
lems with Wilson Quarks 

The location of the finite temperature phase transition Kt is identified by a 
sudden change of physical observables such as the plaquette, the Polyakov 
line and screening hadron masses. (A more precise determination of the 
location will be given by the maximum point of the susceptibility of a physical 
quantity such as the Polyakov loop. However, our statistics is not high 
enough for it.) See Figs. |^ and ^ for the case of Np = 2 at /3 = 5.0 and 4.5. 



Our numerical results of Kt are summarized in Table |T2|. Results of Kt for 
Np = 2 at Nt = 4: and 6 obtained by us and other groups |T^, ^ |5D|, ^ 
are compiled in Fig. ^. (Results for Np = 3 will be discussed in Sec. ||.) 

We expect, at least near the continuum limit, that as the quark mass 
increases from the chiral limit, the transition becomes weaker with the quark 
mass and it becomes strong again when the quark mass is heavy enough 
to recover the first order transition of the SU(3) gauge theory. The MILC 
collaboration performed a systematic study of the transition at various K 
and (3 and found that, contrary to the expectation, when we decrease K 
from the chiral limit K^ on an Nt = A lattice, the Kt transition becomes once 



very strong at ii" ~ 0.18 and becomes weaker again at smaller K [M. On a 



lattice with Nt = 6 they even found a first order transition at K = 0.17 

0.19 m. 



Looking at the phase diagram shown in Fig. ^ closely, we note that the 
Kt lines initially deviate from the Kc line and then approach the K^. line at 
(3 ~ 4.8 and K ~ 0.18 for A^"* = 4 and at /3 ~ 4.8 — 5.2 and K ~ 0.17 — 
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0.19 for Nt = 6, contrary to the naive expectation that they monotonously 
deviate from the Kc line. The points where strong transitions occur are 
just in the region where the Kt hues approach the K^. line. Therefore, it is 
plausible that the strong transition at intermediate values of is a result of 
lattice artifacts caused by this unusual relation of the Kf and K^. lines [^. 
This unusual relation is probably due to the sharp bend of the Kc line at 
/3 ~ 5.0 which is caused by the cross-over phenomenon between weak and 
strong coupling regions of QCD. Our recent study indeed shows that, with 
an improved lattice action, the distance between the the Kc and Kt lines 
becomes monotonically large when we decrease K and, correspondingly, the 
Kt transition becomes rapidly weak as we decrease K from the chiral limit 
[ P^ . Also the unexpected Nt dependence of in the deconfining phase at 
small P discussed in the previous section, is removed with the same improved 
lattice action. 

The appearance of the lattice artifacts implies that we have to be cautious 
when we try to derive the conclusions in the continuum limit from the numer- 
ical results at finite (3. We also note that Nt = 4 is far from the continuum 
limit and therefore we should take with reservation, in particular, quantita- 
tive values in physical units which are quoted in the following. We, however, 
note that the PCAC relation oc expected from chiral symmetry in 
the confining phase is well satisfied even in the strong coupling region and 
therefore we expect that qualitative feature of the chiral transition such as 
the order of the transition does not affected by lattice artifacts. We certainly 
have to check in future that the conclusions in this article are also satisfied 
when an improved action is adopted. 



7 Numerical Results for Chiral Transitions 



As discussed in Sec. P75| , the chiral transition can be studied along the Kc line 
at the crossing point of the Kt and Kc lines, which we denote as the chiral 
transition point jSct- We first address ourselves to the problem of whether the 
chiral limit of the finite temperature transition exists at all. We then study 
the order of the chiral transition. 



In a previous paper ||T2[ we showed that, when Np > 7, there is a bulk 
first order phase transition at /5 = which separates the confining phase at 
small K from a deconfining phase near the chiral limit at K = 1/4. This 
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implies that the Kt hne does not cross the Kc hne at finite /3 for any Nt. On 
the other hand, when Np < 6, the chiral hmit belongs to the confining phase 
at /5 = 0, which implies that there is a crossing point somewhere at finite (3 
for the case Np < 6. 



7.1 On-Kc method 



In order to identify the crossing point (3ct and study the order of the chiral 
transition there, we take the strategy of performing simulations on the Kc 
line starting from a value of f3 in the deconfining phase and reducing (3. We 
call this method "on- i^c" simulation method. The number of iterations A^inv 
needed for the quark matrix inversion, in general, provides a good indicator 



to discriminate the deconfining phase from the confining phase [p!4| , |33|| . The 
use of A^inv as an indicator is extremely useful on the Kc line, because A^inv 
is enormously large on the Kc line in the confining phase, while it is of 
order several hundreds in the deconfining phase. Therefore there is a sudden 
drastic change of A^inv across the boundary of the two phases. This difference 
is due to the fact that there are zero modes around Kc in the confining phase, 



while none exists in the deconfining phase |T^, Q : We have checked this 
difference for the existence of zero modes in various cases discussed below 
and conclude that the difference of iVinv is not a numerical artifact. 

In the deconfining phase on the Kc line, we measure physical observables 
such as the Polyakov loop, the plaquette and hadron screening masses, as 
usual, after thermalization. From the behavior of physical quantities toward 
f3ct, we are able to study the nature of the chiral transition. In the confining 
phase, on the other hand, it is hard to make the system on the Kc line ther- 
malized due to the enormously large A^inv we encounter in the configuration 
generation. In this case, we only obtain at most bounds for several physi- 
cal quantities by measuring the molecular dynamic time evolution of them 
starting a hot state or a mix state. Although it is unsatisfactory that we can- 
not obtain expectation values for physical quantities in the confining phase, 
the on-i^c method is very powerful to identify the critical point because the 
difference between the two phases is clear already with short time-histories. 
We also check that the crossing point thus determined is consistent with a 
linear extrapolation of the line Kf toward the chiral limit. 
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7.2 Chiral transition for A^^^ = 2 

For the case of QCD with two flavors, studies of an effective a model |1|, 
imply that the order of the chiral transition depends on the the strength of 
the Ua(1) anomaly term at the transition temperature. When the strength is 
zero, it is of first order. However, if the strength of the anomaly term in the 
effective a model is non-zero at the starting point of renormalization transfor- 
mation, it is likely that the effective action is attracted to a 0(4) symmetric 
fixed point under renormalization group transformation . Therefore, it is 
plausible that the chiral transition is of second order. 

Our main results of the measurements for Np = 2 are summarized in 



Tables 13 — 15. 



Let us first discuss the results at Nt = 4. In order to confirm the existence 
of the crossing point, we take the largest (farthest) values of Kc for on- 
Kc simulations, that is, Kc{m'^) for Np = 2 in Table ^ and interpolated 
ones. As discussed previously, Kc{m\) in general depend on the value of 
Nt- However, the differences between those on the iVj = 4 and 8 lattices are 
within numerical uncertainties as shown Table Therefore, we take the 
stringent condition to verify the existence of the crossing point, taking the 
farthest values of K^. 

When we take into account the structure of Kcijn^] T 7^ 0) that it sharply 
turns back at finite /?, we may hit the upper part of it by taking the largest 
values of for the "on-iiTc" method. This, however, does not affect the 
conclusion that the Kt line crosses the Kc line. Our estimates for the value 
of (3ct in this case will be slightly underestimated (cf. Fig. 0). This comment 
applies also for Np = 3 and 6. 

We first perform on-i^'c simulations by the R algorithm to identify the 
crossing point, because it is very time consuming to perform simulations 
with the HMC algorithm due to a low acceptance rate on the Kc line in 
the confining phase. We find that when /5 > 4.0, A^inv stays around several 
hundreds, while for (3 < 3.9 it increases with r and exceeds several thousands 
(see Fig. ^) and in accord with this behavior the plaquette, the Polyakov loop 
and m^r decrease rapidly toward those in the confining phase. Therefore we 
identify the crossing point at /5cf ~ 3.9 — 4.0. This /3ct is consistent with a 
linear extrapolation of the Kf line as is shown in Fig. |^. 

Then we repeat on-Kc simulations by the HMC algorithm for /? > 4.0 in 
order to measure physical observables. The time histories for A^inv at /3 > 4.0 
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plotted in Fig. |^ are obtained with the HMC algorithm, which are similar to 
those with the R algorithm. The Ar should be taken small near j3ct in order 
to keep the acceptance rate reasonably high (for (3 = 4.0, 4.1 and 4.2 we 
use Ar = 0.002, 0.005 and 0.005 to get acceptance rates 0.91, 0.79 and 0.93, 
respectively). The value of thus obtained decreases smoothly toward 
zero as the chiral transition is approached and is consistent with zero at the 
estimated Pet (see Fig. 

We find no two-state signals around f3ct- This is in sharp contrast with the 
Np = 3 and 6 cases where we find clear two-state signals at Pet, as discussed 
below. This, together with the vanishing toward Pd, indicates that the 
chiral transition is continuous (second order or crossover) for Np = 2. 

The results from on-Kc simulations on the Nt = 6 lattice are similar to 
those on the Nt = 4 lattice. The estimated transition point is Pd ~ 4.0 — 
4.2. The value of listed in Table |T1| and plotted in Fig. 0, again decreases 
toward zero as P approaches Pet- For A^^ = 18 with the spatial size 18^ x 24, 
we previously found that the transition is at Pd ~ 4.5 — 5.0 [§]. Although 
the spatial size is not large enough, this result suggests that the shift of Pd 
with Nt is very slow. 



7.3 Chiral transition for A^^ = 3 

Main results of measurements for Np = 3 are summarized in Tables |l^ and 
P!7| . The phase diagram for Np = 3 obtained from our simulations at /3=4.0, 
4.5, 4.7, 5.0 and 5.5 is shown in Fig. p. We find that the Kt line linearly 
approaches to the Kc line. In order to confirm the existence of the crossing 
point by on-Kc simulations, we take the largest (farthest) that is Kc{ml) 
for Np = 2 at /9's we have studied, since this is the most stringent condition 
for the existence of Pd- We use them and interpolated values for on-i^'c 
simulations here. For Np = 6 discussed in the next subsection, we interpolate 
these values of Kc with Kc = 0.25 at P = 0. Note that the differences of KcS 
for Np = 2, 3 and 6 are of the same magnitude of numerical uncertainties of 
Kc. 

Fig. 1^ shows A^inv as a function of the molecular-dynamics time r for 
several values of /3's. When P > 3.1, A^in^ is of order of several hundreds, 
while when P < 2.9, A^inv shows a rapid increase with r. At /9 = 3.0 we see 
a clear two-state signal depending on the initial condition: For a hot start, 
A^inv is quite stable around ~ 800 and is large (~ 1.0). On the other 
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hand, for a mix start, A^inv shows a rapid increase with r and exceeds 2,000 
in r ~ 20, and in accordance with this, decreases with r. 

The value of is plotted in Fig. [TO. At /? = 3.0 we have two values 



for depending on the initial configuration. The larger one obtained for 
the hot start is of order 1.0, which is a smooth extrapolation of the values at 
/? ~ 3.1 - 3.2. The smaller one is an upper bound for for the mix start. 

We note that the result of Pet ~ 3.0 is consistent with an extrapolation 
of Kt points listed in Table |1^ as is shown in Fig. |^. (The nature of the 
transition Kf off the chiral limit is discussed in Sec. p.) Thus we identify the 
crossing point at Pet ~ 3.0(1). With the clear two-state signal we conclude 
that the chiral transition is of first order for Np = 3. 

7.4 Chiral transition for A^^^ = 6 



Our previous study at /3 = |T2| shows that for Np = 7 there is no crossing 
point of the Kc and Kt lines and that A^^ = 6 is the largest number of flavors 
for which a crossing point exists. Main results of measurements for Np = 6 



are summarized in Table IS. Overall features of the transition obtained from 



numerical simulations for Np = 6 are very similar to those for Np = 3 except 



for the location of P^, which moves to a smaller P as expected. Fig. [11] shows 
that A^inv on the Kc line stays at several hundreds for P > 0.4 and for a hot 
start at P = 0.3. On the other hand, iVjnv grows rapidly with r and exceeds 
5,000 for P < 0.2 and for a mix start aX P = 0.3. In accord with this, we have 
two values of at /? = 0.3 (cf. Fig. |12D. Therefore we identify the crossing 



point at Pet ~ 0.3(1) and conclude that the chiral transition is of first order 
for Np = 6. This Pd is consistent with a linear extrapolation of the Kt line 
(cf. Table [ID. 

For QCD with Np > 3, Pisarski and Wilczek predicted a first order chiral 
transition from a renormalization group study of an effective a model [^. Our 
results for Np = 3 and 6 are consistent with their prediction. 



8 Influence of the Strange Quark 

In the previous section, we have seen that the chiral transition is consistent 
with a second order transition for Np = 2, while it is of first order for Np > 3, 
both in accordance with theoretical expectations. Off the chiral limit, we ex- 
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pect that the first order transition for Np > 3 smoothens into a crossover 
at sufficiently large rrig. In this way the nature of the transition sensitively 
depends on Np and rrig. Therefore, in order to study the nature of the tran- 
sition in the real world, we should include the strange quark properly whose 
mass rUs is of the same order of magnitude as the transition temperature 
~ 100 — 200 MeV. 

In a numerical study we are able to vary the mass of the strange quark. 
When the mass of the strange quark is reduced from infinity to zero with up 
and down quarks fixed to the chiral limit, the nature of the transition must 
change from continuous to first order at some quark mass m*. Assuming 
that the chiral transition is of second order for Np = 2 (i.e. = oo), this 
point at m* is a tricritical point 0. The crucial question is whether the 
physical strange quark mass is larger or smaller than m*. Studies with an 
effective linear a model suggest a crossover for the case of realistic quark 
masses in meanfield approximation and in a large 1/Np approximation [^G], 



37| , while the possibility of a weakly first order transition is not excluded 



when numerical errors in the calculation of basic parameters are taken into 



account 37 



8.1 Nf = 3 

Let us first discuss the case of the degenerate Np = 3: Ku = Kd = Kg = K. 
As we have already discussed the chiral transition previously, we are mainly 
interested in the transition for the massive quarks. In order to find the 
transition points we perform simulations at /3=4.0, 4.5, 4.7, 5.0 and 5.5. The 



results for physical quantities are plotted in Figs. |I3] — |T^. The transition 
points identified by a sudden change of physical observables are given in 
Table and plotted in Fig. |. We note that the Kt line for A^^ = 3 at A^^t = 4 
locates sufficiently far from the points where the Kc line bends rapidly. This 
situation is quite different from the Np = 2 case where the unusual relation 
between the Kt line and line causes the lattice artifacts. Therefore, we 
expect that these lattice artifacts are small in the Np = 3 case. 

In the previous section we have seen that the transition is of first order 
in the chiral limit Kc = 0.235 at f3 = 3.0 for Np = 3. For phenomenological 
applications, it is important to estimate the critical value of the quark mass 
m'^"^ up to which the first order phase transition persists. 

We observe clear two state signals at /3 = 4.0, 4.5 and 4.7, while for /3 = 5.0 
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and 5.5 no such signals have been seen: The simulation time history of the 
plaquette at /3 = 4.7 on a 12^ x 4 lattice is plotted in Fig. |l8|(a). The confining 
and deconfining phases coexist over 1,000 trajectories at K = 0.1795 and, in 
accordance with this, we find two-state signals also in other observables such 
as the plaquette and the pion screening mass (cf. Fig. |l^). From them 
we conclude that the transition at i^' = 0.1795(5) and (3 = 4.7 is first order. 
On the other hand, the time history of the plaquette at /3 = 5.0 shown in 
Fig. |l5(b) suggests that the transition is a crossover there. 

At the transition point (in the confining phase) of f3 = 4.7 the value of 
ruqa is 0.175(2) and mT^/nip = 0.873(6). The results of the hadron spectrum 
in the range of /3 = 3.0 - 4.7 for Np = 2 and 3 (cf. Fig. |l^) indicate that 
the inverse lattice spacing estimated from the rho meson mass is almost 
independent on /3 in this range and a^^ ~ 0.8 GeV. (Hereafter we use 
determined from rUp in the chiral limit.) Therefore we obtain a bound on the 
critical quark mass m^^** > 140 MeV, or equivalently (m^/mp)^''** > 0.873(6). 
It should be noted that the physical strange quark mass determined from 

= 1020 MeV, using the data shown in Fig. turns out to be rris ~ 150 
MeV in this (3 range with our definition of the quark mass. 

We note that these values for the critical quark mass are much larger 
than those with staggered quarks where rrf^'^^a = 0.025 — 0.075 [^, ^ 
(m^"* -10 — 40 MeV using a"^ ~ 0.5 GeV at /3 = 5.2 hr Np = 2 ^) 
which means that {m.,^/mpY^^^ ~ 0.42 — 0.58 (using the results of meson 
masses for Np = 4 at /3 = 5.2 [^, because the data for Np = 3 are not 
available). 



8.2 Nf = 2 + 1 

Now let us discuss a more realistic case of massless up and down quarks 
and a light strange quark [Np = 2 + 1). Main results of measurements are 
summarized in Tables ^ — |2T]. Our strategy to study the phase structure is 
similar to that applied in Sec. |^ for the investigation of the chiral transition 
in the degenerate quark mass cases, which we called the on- i^c method. We 
set the value of masses for the up and down quarks niud to zero {Kud = Kc) 
and fix the strange quark mass rris to some value, and make simulations 
starting from a value of (3 in the deconfining phase and reducing the value of 
p. When u and d quarks are massless, the number of iteration Ni^v needed 
for the quark matrix inversion (for u and d quarks) is enormously large in 
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the confining phase, while it is of order of several hundreds in the deconfining 
phase. The values which we take for are given in Table 22. They are the 
vanishing point of extrapolated for Np = 2 and interpolated ones. We 
have used those for Np = 2, because we have the data most in this case, 
and the difference between that for Np = 2 and 3 is of the same order of 
magnitude as the difference due to the definition of Kc (cf. discussions in 
Sec. 0). 

We study two cases of ~ 150 MeV and 400 MeV. From the value of 
~ 0.8 GeV and an empirical rule m^a ~ (2/3)(l/i^ — 1/Kc) satisfied for 
Np = 2 and 3 in the /5 region we have studied (cf. Fig. |T9|), we get the values 



for Kc: shown in Table 22 



In order to confirm that our choice of parameters for the case ~ 
150 MeV is really close to the physical values, we have also made a zero- 
temperature spectroscopy calculation for the Np = 2 + 1 case at /5 = 3.5 
on an 8^ x 10 lattice. Keeping Kg = 0.2017 (m^ ~ 150 MeV), we vary Ky^^ 
from 0.195 to 0.210 in steps of 0.005. Taking the chiral limit of Kud, we 
obtain = 903(38) MeV from the rho meson mass [rripa = 0.853(36) at 
Kc = 0.2227, where Kc is determined by a linear extrapolation of m^a^ in 
terms of 1/K). The mass of the meson at the simulation point turns out to 
be 1.03(5) GeV which should be compared with the physical value 1.02 GeV. 
Thus the hopping parameter chosen for ~ 150 MeV corresponds to the 
physical strange quark mass, in this sense. As far as we consider the meson 
sector the numerical results for mass ratio do not differ so much from the 
physical values. However, we emphasize one caveat here. The nucleon-rho 
mass ratio /nip turns out to be 2.0(1) which is the same as the result 2.0 
in the strong coupling limit and is much larger than the physical value 1.22. 
This implies that /5=3.5 is far from the continuum limit. 

The simulation time history of A^inv on the 8^ x 10 spatial lattice is plotted 
in Fig. ^OKa) for the case of 150 MeV. When [3 > 3.6, A^inv is of order of several 
hundreds, while when < 3.4, iYjnv shows a rapid increase with r. At /3 = 3.5 
we see a clear two-state signal depending on the initial condition: For a hot 
start, iVinv is quite stable around 900 and m\ is large (~ 1.0 in lattice units). 
On the other hand, for a mix start, A^inv shows a rapid increase with r and 
exceeds 2,500 in r ~ 10, and in accord with this, the plaquette and m\ 
decreases with r as shown in Fig. ^(b) for the plaquette. For the case of 
400 MeV a similar clear two-state signal is observed at /? = 3.9 both on the 
8^ X 10 and 12^ spatial lattices (cf. Fig. [2T|). The values of versus (3 are 
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plotted in Fig. together with those in the case of degenerate Np = 3 on 
the Kc line. At (3 = 3.5 for the case of 150 MeV and at /? = 3.9 for the case of 
400 MeV, we have two values for depending on the initial configuration. 
The larger ones of order 1.0 are for hot starts, while the smaller ones are 
upper bounds for mix starts. These results imply that m* > 400 MeV in our 
normalization for quark masses. 



Following the Columbia group |39|, we summarize our results about the 



nature of the QCD transition at A'f = 4 as a function of rriud and nis in Fig. 23 



together with theoretical expectations [jT], |], |4^ assuming that the chiral 
transition is of second order for Np = 2. Clearly the point which corresponds 
to the physical values of the up, down and strange quark masses measured 
by m^/mp and m.,t/mp exists in the range of the first order transition. If 
this situation persists in the continuum limit, the transition for the physical 
quark masses is of first order. 

The Columbia group studied the influence of the strange quark for the 
case of staggered quarks |^. Their result shows that no transition occurs 



at m^a = m^a = 0.025, nisa = 0.1 (m^ = ~ 12 MeV, ~ 50 MeV 
using ~ 0.5 GeV). Their zero-temperature values for mK/frip and niT^/mp 
obtained at this simulation point suggest that this value for iris is smaller 
than its physical value and those for ruu and irid are larger than their physical 
values. This implies that the transition in the real world is also a crossover, 
unless the second order transition line, which has a sharp iriud dependence 



near m* as shown in Fig. |2^ crosses between the physical point and the 
simulation point. 

Although both staggered and Wilson simulations give phase structures 
qualitatively consistent with theoretical expectations JI], ^ we note that 
Wilson quarks tend to give larger values for critical quark masses (measured 
by m^/mp etc.) than those with staggered quarks. This leads to the differ- 
ence in the conclusions about the nature of the physical transition. However, 
since the deviation from the continuum limit is large in both of the studies at 
Nt = 4, we certainly should make a calculation with larger Nt |^ or using 
an improved action to get closer to the continuum limit and to obtain 
a definite conclusion about the nature of the QCD transition. With Wilson 
quarks using the standard gauge action, however, A^^ should be enormously 
large (> 18) [|| in order to avoid the lattice artifacts discussed in Sec. ^ Im- 
provement of the lattice action will be essential especially for Wilson quarks. 
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9 Conclusions 



We have studied the nature of finite temperature transitions near the chiral 
hmit for various numbers of fiavors {Np = 2, 3, and 6) and also for the case 
of massless up and down quarks and a hght strange quark {Np = 2 + 1), 
mainly on lattices with Nt = 4, using the Wilson formalism of quarks on the 
lattice. 

We have found that the chiral transition is continuous (second order or 
crossover) for Np = 2, while it is of first order for Np = 3 and 6. These 
results are in accordance with theoretical predictions based on universality 
[0, 0]. Our results with Wilson quarks are also consistent with those with 
staggered quarks ||4^ . 

Our results for QCD with a strange quark as well as up and down quarks 
obtained on Nt = 4 lattices are summarized in Fig. ^ Clearly, the point 
which corresponds to the physical values of the up, down and strange quark 
masses measured by / rrip and rriT^/mp, marked with star in Fig. exists in 
the range of first order transition. If this situation persists in the continuum 
limit, the transition for the physical quark masses is of first order. 

We have found that Wilson quarks tend to give larger values for critical 
quark masses (measured, for example, by m^/rrip and m.„/mp ) than those 
with staggered quarks. This leads to the difference in the conclusions about 
the nature of the physical transition. Because the deviation from the con- 
tinuum limit is large on the A'^t = 4 lattices, we certainly should make a 
calculation with larger A^^ or with an improved action |^ in order to get 
closer to the continuum limit and to obtain a definite conclusion about the 
nature of the physical QCD transition, by resolving the discrepancy between 
Wilson and staggered quarks for the conclusions. Studies with an improved 
gauge action and the Wilson quark action are in progress. 
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(3 




Ar 




UllL^J. Ill 


alffo. 

o 


iVinv 


phase 





0, 


.2 


0.02 


1132 


500 


H-CR 


37 


c 





0, 


.21 


0.01 


1005 


500 


H-CR 


48 


c 





0, 


.22 


0.01 


1041 


500 


H-CR 


45 


c 





0, 


.23 


0.01 


700 


500 


H-CR 


95 


c 


3 


0, 


.18 


0.01 


250 


100 


H-CR 


37 


c 


3 


0, 


.19 


0.01 


150 


100 


H-CR 


35 


c 


3 





.2 


0.01 


160 


100 


H-CR 


48 


c 


3.5 


0, 


.175 


0.01 


160 


100 


H-CR 


27 


c 


3.5 





.185 


0.01 


160 


100 


H-CR 


34 


c 


3.5 


0, 


.195 


0.01 


160 


100 


H-CR 


46 


c 


4 





.17 


0.02 


1650 


500 


H-CR 


15 


c 


4 


0, 


.18 


0.02 


2188 


1000 


H-CR 


18 


c 


4 


0, 


.19 


0.02 


1550 


500 


H-CR 


23 


c 


4 


0, 


.2226 


0.002 


50 


24 


H-CG 


1054 


d 


4.1 


0, 


.2211 


0.005 


92 


50 


H-CG 


781 


d 


4.2 


0, 


.2195 


0.005 


206 


100 


H-CG 


430 


d 



Table 1: Table of job parameters for Np = 2 simulations performed on an 
8^ X 10 X 4 lattice. Data marked with f are taken from our previous simulation 
[|l^ performed on an 8^ x 20 x 4 lattice. The column "algo." is for the 
algorithm used for update (HMC or R) and for quark matrix inversion (CR 
or CG). A^inv is an average number of iterations needed for the quark matrix 
inversion. Errors for A^inv are in most cases about 1%. The last column is for 
the initial and final phases (c: the low temperature confining phase, d: the 
high temperature deconfining phase, and m: mix state), where parentheses 
mean that the system is not completely thermalized. 
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p 


K 


At 


Ttot 


^therm 


algo. 


N- 


phase 


4.3 


0.165 


0.02 


520 


320 


H-CR 


23 


c 


4.3 


0.175 


0.01 


490 


290 


H-CR 


28 


c 


4.3 


0.185 


0.01 


400 


200 


H-CR 


39 


c 


4.3 


0.205 


0.008 


460 


250 


H-CR 


250 


d^c 


4.3 


0.207 


0.005 


16 




H-CG 




c 


4.3 


0.207 


0.005 


30 




H-CG 




d^(c) 


4.3 


0.208 


0.005 


38 




H-CG 




c 


4.3 


0.208 


0.005 


45 




H-CG 




d^(c) 


4.3 


0.21 


0.005 


150 


50 


H-CG 


820 


d 


4.3 


0.218 


0.01 


196 


100 


H-CG 


338 


d 


4.5 


0.16 


0.02 


500 


300 


H-CR 


25 


c 


4.5 


0.17 


0.01 


580 


300 


H-CR 


34 


c 


4.5 


0.18 


0.01 


530 


300 


H-CR 


42 


c 


4.5 


0.195 


0.01 


310 


100 


H-CR 


92 


c 


4.5 


0.2 


0.005 


175 


135 


H-CR 


280 


c 


4.5 


0.202 


0.008 


700 


300 


H-CG 


473 


d 


4.5 


0.205 


0.01 


190 


100 


H-CG 


314 


d 


4.5 


0.2143 


0.01 


197 


100 


H-CG 


209 


d 


5 


0.14 


0.02 


500 


300 


H-CR 


17 


c 


5 


0.15 


0.02 


520 


300 


H-CR 


20 


c 


5 


0.16 


0.02 


600 


300 


H-CR 


24 


c 


5 


0.17 


0.01 


540 


300 


H-CR 


41 


d^c 


5 


0.18 


0.01 


640 


200 


H-CG 


169 


c^d 


5 


0.19 


0.01 


720 


300 


H-CG 


132 


d 


5 


0.1982 


0.01 


761 


300 


H-CG 


118 


c^d 



Table 1: Continued. 
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p 


K 




'^tot 


^therm 


algo. 


AT 

i * inv 


phase 


5.25 


0. 


1 


0.01 


r oo 

520 


300 


TT /^T> 

H-CR 


12 


c 


5.25 


0. 


11 


0.01 


600 


300 


TT /^T^ 

H-CR 


13 


c 


r or 

5.25 


0. 


12 


r\ ^\-^ 

U.Ul 


600 


300 




15 


c 


5.25 


0. 


13 


r\ r\-\ 

0.01 


r" r\ 

560 


300 


TT / ^ T > 

H-CR 


1 r' 

17 


c 


5.25 


0. 


,14 


0.01 


r o 

580 


300 


TT /'I T"> 

H-CR 


20 


c 


5.25 


u. 


,io 


0.01 


520 


300 


TT /^T^ 

H-CR 


o r 

25 


c 


o.zo 


0. 


155 


U.Ul 


ozu 


OUU 


W PR 

n- 




a — >L 


5.25 


0. 


16 


0.01 


540 


300 


H-CR 


39 




5.25 


0. 


165 


0.01 


600 


300 


H-CG 


121 


d 


5.25 


0. 


175 


0.01 


610 


300 


H-CG 


118 


d 


5.25 


0. 


18 


0.01 


640 


300 


H-CG 


111 


d 


5.5t 


0. 


15 


0.025 


2500 


800 


H-CR 


17 


d 


5.5t 


0. 


16 


0.025 


1572 


500 


H-CR 


37 


d 


5.5t 


0. 


1615 


0.025 


1532 


500 


H-CR 


43 


d 


5.5t 


0. 


163 


0.025 


1458 


500 


H-CR 


53 


d 


6 


0. 


15 


0.01 


427 


200 


H-CG 


73 


d 


6 


0. 


1524 


0.01 


230 


150 


H-CG 


78 


d 


6 


0. 


155 


0.01 


427 


200 


H-CG 


80 


d 


6 


0. 


16 


0.01 


400 


200 


H-CG 


83 


d 


10 


0. 


13 


0.01 


351 


200 


H-CG 


48 


d 


10 


0. 


14 


0.01 


400 


200 


H-CG 


78 


d 


10 


0. 


15 


0.01 


338 


200 


H-CG 


60 


d 



Table 1: Continued. 
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/3 


K 


Ar 


Ttot 


'''therm 


algo. 


iVinv 


phase 


4.2 


0.2195 


0.00125 


56 


30 


H-CG 


1119 


d 


4.3 


0.2183 


0.002 


138 


40 


H-CG 


863 


d 


4.4 


0.2163 


0.005 


160 


30 


H-CG 


678 


d 


4.5 


0.2143 


0.008 


130 


80 


H-CG 


505 


d 


5 


0.1982 


0.01 


224 


100 


H-CG 


160 


d 


5.02 


0.16 


0.01 


560 


300 


H-CR 


27 


c 


5.02 


0.17 


0.01 


560 


300 


H-CR 


36 


c 


5.02 


0.18 


0.01 


180 


100 


H-CG 


143 


c 


5.02 


0.18 


0.01 


210 


100 


H-CG 


529 


d 



Table 2: The same as Table |T| for Np = 2 simulations performed on a 12'^ x 6 
lattice. 



p 


K 


Ar 


Ttot 


''therm 


algo. 




phase 


4.5 


0.16 


0.02 


500 


300 


H-CR 


23 


c 


4.5 


0.17 


0.01 


540 


300 


H-CR 


29 


c 


4.5 


0.18 


0.01 


540 


300 


H-CR 


35 


c 


5.5t 


0.15 


0.025 


2050 


1000 


H-CR 


8 


c 


5.5t 


0.155 


0.02 


1600 


500 


H-CR 


23 


c 


6 


0.1524 


0.01 


230 


150 


H-CG 


78 


d 



Table 3: The same as Table |l| for Np = 2 on an 8^ x 10 lattice. Data marked 
with t are taken from Ref. |]TE[ obtained on an 8^ x 20 lattice. 
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K 


Ar 


Ttot 


^therm 


algo. 




phase 


2.5 


0.2381 


0, 


.01 


8 




R-CG 


-2300 


d-(c) 


2.7 


0.2369 


0, 


.01 


10 




R-CG 


-2300 


d-(c) 


2.8 


0.2364 


0, 


.01 


12 




R-CG 


>1900 


d-(c) 


2.9 


0.2358 


0, 


.01 


28 




R-CG 


-2300 


d-(c) 


3 


0.205 


0, 


.01 


280 


170 


R-CR 


64 


c 


3 


0.205 


0, 


.01 


202 


100 


R-CR 


64 


c 


3 


0.215 


0, 


.01 


190 


100 


R-CR 


117 


c 


3 


0.225 


0, 


.005 


75 




R-CR 


-563 


d-(c) 


3 


0.23 


0, 


.0025 


18 




R-CG 




d-(c) 


3 


0.2352 


0, 


.01 


23 




R-CG 


-2300 


m^(c) 


3 


0.2352 


0, 


.01 


68 




R-CG 


-2300 


c-^(c) 


3 


0.2352 


0, 


.01 


159 


100 


R-CG 


851 


d 


3.1 


0.2341 


0, 


.01 


160 


50 


R-CG 


650 


d 


3.2 


0.2329 


0, 


.01 


114 


50 


R-CG 


556 


d 


3.2 


0.2329 


0, 


.01 


169 


100 


R-CG 


504 


d 


4 


0.18 


0, 


.01 


520 


300 


R-CR 


35 


c 


4 


0.19 


0, 


.01 


520 


300 


R-CR 


47 


c 


4 


0.2 


0, 


.01 


391 


200 


R-CR 


84 


d^c 


4 


0.205 


0, 


.01 


320 


200 


R-CG 


351 


d 


4 


0.21 


0, 


.01 


308 


200 


R-CG 


247 


d 


4 


0.2226 


0, 


.01 


320 


200 


R-CG 


188 


d 


4.5 


0.16 


0, 


.01 


500 


300 


R-CR 


25 


c 


4.5 


0.17 


0, 


.01 


542 


300 


R-CR 


30 


c 


4.5 


0.18 


0, 


.01 


545 


300 


R-CR 


40 


d^c 


4.5 


0.185 


0, 


.01 


534 


300 


R-CR 


51 


d^c 


4.5 


0.186 


0, 


.01 


301 


150 


R-CR 


56 


c 


4.5 


0.1875 


0, 


.01 


191 


100 


R-CR 


82 


c 


4.5 


0.1875 


0, 


.01 


181 


100 


R-CG 


248 


d 


4.5 


0.189 


0, 


.01 


207 


100 


R-CG 


214 


d 



Table 4: The same as Table I for iVi^ = 3 on an 8^ x 10 x 4 lattice. 
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K 


Ar 
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^ * mv 


phase 


4.5 


0.19 


0.01 


336 


200 


R-CG 


200 


d 


4.5 


0.2 


0.01 


394 


200 


R-CG 


158 


d 


4.5 


0.205 


0.01 


190 




R-CG 


142 


d 


4.5 


0.2143 


0.01 


101 




R-CG 


132 


d 


5 


0.13 


0.01 


313 


150 


R-CR 


49 


c 


5 


0.14 


0.01 


275 


150 


R-CR 


20 


c 


5 


0.15 


0.01 


310 


150 


R-CR 


22 


c 


5 


0.16 


0.01 


324 


150 


R-CR 


27 


c 


5 


0.165 


0.01 


373 


150 


R-CR 


65 


c 


5 


0.165 


0.01 


202 


150 


R-CG 


48 


d— s>c 


5 


0.166 


0.01 


120 




R-CG 




d^(c) 


5 


0.166 


0.01 


264 


150 


R-CR 


35 


c 


5 


0.167 


0.01 


145 




R-CR 




c^(d) 


5 


0.167 


0.01 


187 


100 


R-CG 


155 


d 


5 


0.17 


0.01 


291 


150 


R-CG 


119 


d 


5.5 


0.1 


0.01 


652 


150 


R-CR 


13 


c 


5.5 


0.11 


0.01 


505 


150 


R-CR 


15 


c 


5.5 


0.12 


0.01 


571 


250 


R-CR 


16 


c 


5.5 


0.125 


0.01 


695 


250 


R-CR 


17 


c— >m 


5.5 


0.1275 


0.01 


676 


100 


R-CR 


18 


d 


5.5 


0.13 


0.01 


364 


150 


R-CR 


18 


c^d 


5.5 


0.135 


0.01 


174 


100 


R-CR 


20 


d^d 


5.5 


0.14 


0.01 


296 


100 


R-CR 


23 


d 


6 


0.08 


0.01 


355 


100 


R-CG 


23 


d 


6 


0.09 


0.01 


194 


100 


R-CG 


27 


d 


6 


0.1 


0.01 


320 


100 


R-CG 


33 


d 


6 


0.11 


0.01 


270 


100 


R-CG 


41 


d 


6 


0.12 


0.01 


244 


100 


R-CG 


51 


d 


6 


0.135 


0.01 


180 


100 


R-CG 


72 


d 



Table 4: Continued. 
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n ni 
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onQ 
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1 nn 
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R PP 
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Cl 


4 5 


n 18 


01 


282 
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R-CR 


40 




4.5 


0.186 


0.01 


230 


100 


R-CR 


56 


c 


4.5 


0.1875 


0.01 


1040 


369 


R-CR 


74 


c 


4.5 


0.1875 


0.01 


1072 


100 


R-CG 


264 


d 


4.5 


0.189 


0.01 


183 


100 


R-CG 


230 


d 


4.5 


0.19 


0.01 


196 


100 


R-CG 


219 


d 


4.7 


0.17 


0.01 


307 


100 


R-CR 


32 


c 


4.7 


0.175 


0.01 


225 


100 


R-CR 


38 


c 


4.7 


0.178 


0.01 


232 


100 


R-CG 


117 


d^c 


4.7 


0.179 


0.01 


335 


100 


R-CR 


48 




4.7 


0.179 


0.01 


253 




R-CG 




d-(c) 


4.7 


0.1795 


0.01 


1035 


100 


R-CR 


50 


c 


4.7 


0.1795 


0.01 


1073 


100 


R-CG 


236 


d 


4.7 


0.18 


0.01 


299 


100 


R-CG 


228 


d 


4.7 


0.18 


0.01 


410 




R-CG 




c-(d) 



Table 5: The same as Table for Np = 3 on a 12^ x 4 lattice. 
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K 


Ar 


''tot 


''"therm 


algo. 


N- 


phase 


5 


0.165 


0.01 


203 


100 


R-CR 


33 


c 


5 


0.166 


0.01 


574 


200 


R-CR 


35 




5 


0.166 


0.01 


405 




R-CG/CR 




d-(c) 


5 


0.16625 


0.01 


570 


200 


R-CR 


37 


c— >m 


5 


0.16625 


0.01 


545 


200 


R-CR 


47 


d— >m 


5 


0.1665 


0.01 


502 




R-CR 




c-(d) 


5 


0.1665 


0.01 


611 


200 


R-CR 


75 


d 


5 


0.167 


0.01 


475 


250 


R-CR 


53 


d 


5 


0.168 


0.01 


419 


100 


R-CR 


104 


d 


5 


0.169 


0.01 


164 


100 


R-CG 


163 


d 


5 


0.17 


0.01 


231 


100 


R-CG 


166 


d 



Table 5: Continued. 



33 



P 


j\ 




Ttot 




algo. 


AT 


pildbc 


n 
u 


u./ 


U.Ui 


OZ 


on 
zu 


R PR 
rL-oix 


'ifi 
Oo 


C 


n 
u 


U.Zi 


n m 

U.Ui 


oz 


on 
zu 




AQ 
4y 


C 


n 
u 


u./z 


U.Ui 


OO 


1 fi 

io 


R PR 


( 


C 


n 
u 


U.ZoO 


n ni 

U.Ui 


A n 
4U 


on 
zu 


R PR 


iOO 


C 


n 1 

U.i 


u.z4yo 


n ni 
U.Ui 


i i 




R PP 


^i^nnn 
>OUUU 




U./ 


n o/io 
u.z4y 


n ni 
U.Ui 


ii 




R PP 


\ cnnn 
>OUUU 




n o 
U.z 


n 0/10QR 

U.z4yoo 


n ni 
U.Ui 


Zo 




R PP 


\ 1 nnnn 
>iUUUU 




n Q 
U.o 


U.z4oO 


n ni 
U.Ui 


iD 




R PP 


\ cnnn 
>OUUU 


m^(c) 


n Q 
U.o 


U.z4oO 


n ni 
U.Ui 


O 

y 




R PP 


\ cnnn 
>OUUU 


m^(c) 


u.o 


U.Z4oO 


n ni 
U.Ui 


Z ( 




R PP 


Rnn 
ouu 


Q 


U.o 


n o/in 


n ni 
U.Ui 


1 R 

io 






\ cnnn 
>OUUU 


m^(c) 


n A 

U.4 


n o/ifi 

U./4o 


n ni 

U.Ui 


on 
zu 


1 n 

iU 


R PP 


c;nn 
ouu 




u.o 


U.zo 


n ni 

U.Ui 







R PP 






u.o 


U.ZoO 


n ni 

U.Ui 


D 




R PP 






u.o 


n OA 

U.Z4 


n ni 

U.Ui 







R PP 




J .( „\ 


u.o 


U.Z40 


n ni 

U.Ui 


Oo 




R PP 


^ ,1 /I nn 

~i4UU 


Q — ^-C 


u.o 


U.z4 / 


n ni 

U.Ui 


OC^ 

zo 


1 K 
iO 


R PP 


440 


(1 


1 


u./ 
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01 
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50 


R-CR 


60 


n 


1 


0.22 


0.01 


115 


55 


R-CR 


80 


c 


1 


0.225 


0.01 


267 


100 


R-CR 


126 


c 


1 


0.23 


0.01 


293 


100 


R-CR 


192 


c 


1 


0.235 


0.01 


40 




R-CG 




d^c 


1 


0.235 


0.005 


112 


60 


R-CG 


970 




1 


0.235 


0.005 


19 




R-CG 




d-(c) 


1 


0.237 


0.005 


42 




R-CG 




d 


1 


0.237 


0.005 


49 




R-CG 




c-(d) 


1 


0.238 


0.005 


28 




R-CG 


440 


d 


1 


0.24 


0.005 


108 


40 


R-CG 


325 


d 


1 


0.245 


0.01 


114 


60 


R-CG 


306 


d 



Table 6: The same as Table |T] for Np = 6 on an 8^ x 10 x 4 lattice. 
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K 


Ar 


Ttot 
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algo. 


A^inv 


phase 


2 


0.24 


0.01 


18 




R-CG 


162 


d 


4 


0.22 


0.01 


15 




R-CG 


88 


d 


4.5 


0.15 


0.01 


71 


61 


R-CR 


21 


c 


4.5 


0.16 


0.01 


38 


28 


R-CR 


27 


c 


4.5 


0.165 


0.01 


60 


50 


R-CR 


32 




4.5 


0.165 


0.01 


60 




R-CR 




d-(c) 


4.5 


0.166 


0.01 


277 


267 


R-CR 


36 


d^c 


4.5 


0.167 


0.01 


193 


183 


R-CR 


36 


c 


4.5 


0.167 


0.01 


159 


149 


R-CR 


105 


d 


4.5 


0.168 


0.01 


152 




R-CG 




c^(d) 


4.5 


0.17 


0.01 


73 




R-CG 




c^(d) 


4.5 


0.18 


0.01 


41 


31 


R-CG 


115 


c^d 


4.5 


0.19 


0.01 


38 


28 


R-CG 


92 


c^d 


4.5 


0.2143 


0.01 


181 


150 


R-CG 


87 


d 



Table 6: Continued. 
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Kud 


Ks 


''"tot 


^therm 


J\Tud 
inv 


IIIV 


phase 


3.2 


0.2329 


0.2043 


15 




-3000 


-250 


d^(c) 


3.4 


0.2306 


0.2026 


20 




-3000 


-290 


d-(c) 


3.5 


0.2295 


0.2017 


9 




-3000 


-260 


m— >(c) 


3.5 


0.2295 


0.2017 


553 


100 


862 


394 


d 


3.6 


0.2281 


0.2006 


153 


100 


622 


344 


d 


3.7 


0.2267 


0.1692 


20 




-2500 


-100 


d-(c) 


3.8 


0.2254 


0.1684 


47 




-2500 


-93 


d-(c) 


3.9 


0.224 


0.1677 


12 




-2500 


-93 


m— ^(c) 


3.9 


0.224 


0.1677 


760 


100 


797 


135 


d 


4 


0.2226 


0.1669 


159 


100 


521 


137 


d 


4 


0.2226 


0.1964 


167 


100 


235 


201 


d 


4.3 


0.218 


0.1643 


159 


100 


229 


130 


d 


5.5 


0.163 


0.15 


376 


208 


119 


97 


d 



Table 7: The same as Table |1| for Np = 2 + lonan8^xl0x4 lattice. The 
molecular dynamics time step is Ar = 0.01. Simulations are performed with 
the R algorithm for updating configurations and with the CG method for 
quark matrix inversions. 



p 




Ks 


T~tot 


^therm 


mv 


^i^nv phase 


3.9 


0.224 


0.1677 


14 




-3000 


—93 m-^(c) 


3.9 


0.224 


0.1677 


398 


100 


999 


139 d 


4 


0.2226 


0.1669 


396 


100 


636 


141 d 



Table 8: The same as Table ^ for Np = 2 + 1 on an 8^ x 10 x 4 lattice. 
Simulations are performed with the R algorithm for updating configurations 
and with the CG method for quark matrix inversions. 
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/? 


Kud 


Ks 


Ttot 




J\Tud 
inv 


mv 


phase 


3.5 


0.195 


0.2017 


196 


100 


46 


58 


c 


3.5 


0.2 


0.2017 


164 


50 


57 


59 


c 


3.5 


0.205 


0.2017 


166 


50 


74 


61 


c 


3.5 


0.21 


0.2017 


158 


40 


109 


64 


c 



Table 9: The same as Table |^ for Np = 2 + 1 on an 8^ x 10 lattice. Simulations 
are performed with the R algorithm for updating configurations and with the 
CR method for quark matrix inversions. 



p 






3.0 


0.235(1) 


0.230(1) 


3.5 


0.230(1) 


0.226(1) 


4.0 


0.223(1) 


0.218(4) 


4.3 


0.218(1) 


0.214(1) 


4.5 


0.214(1) 


0.210(1) 


6.0 




0.1564(1) 


10.0 




0.1396(1) 



Table 10: The chiral limit Kc for Np = 2 determined on an 8^ x 10 x 4 lattice. 
The results for j3 = 3.0 — 4.5 are determined by = and rriq = 0, where 
values of and in the confining phase are linearly extrapolated in 1/K. 
The results for f3 = 6.0 and 10.0 are determined from an interpolation of mq 
in the deconfining phase. 
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Nt 


= 4 


Nt = S 


Nf 


Kc{ml) 






2 


0.214(1) 


0.210(1) 


0.212(1) 0.209(1) 


3 


0.210(1) 


0.204(1) 




6 


0.205(2) 


0.200(1) 





Table 11: The chiral limit K^. aX, (3 = A.h determined by = and ruq = 0, 
where values of and rriq in the confining phase are linearly extrapolated 
in 1/K using data from = 0.16 — 0.18 for Nf = 2 and 3 and K = 0.15 
— 0.165 for Nf = 6 (because Kt = 0.167(1) for iV^. = 6 at iVi = 4). The 
spatial lattice size is 8^ x 10. 





Nf = 2 




Nf = S 




A^^ = 6 


/? 


Kt 


/? 


Kt 




Kt 


4.3 


0.207-0.210 


3.0 


> 0.230 


0.5 


0.245-0.2475 


4.5 


0.200-0.202 


4.0 


0.200-0.205 


1.0 


0.235-0.237 


5.0 


0.170-0.180 


4.5 


0.186-0.189 


4.5 


0.166-0.168 


5.25 


0.160-0.165 


4.5* 


0.186-0.189 










4.7* 


0.179-0.180 










5.0 


0.166-0.167 










5.0* 


0.166-0.1665 










5.5 


0.125-0.130 







Table 12: Finite temperature transition Kt for Nf = 2, 3 and 6 obtained on 
an 8^ X 10 X 4 lattice (data with * obtained on a 12^ x 4 lattice). For Nf = 2 



at /5 = 5.0, the data by the MILC collaboration ^6| give a more precise value 
0.177 — 0.178 for Kt (cf. Fig. §. 
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K 


plaquette 


rolyak 


ov 






iTTlqa 


VTLpa 





0.2 


0.0088(1) 


0.0367 


.3) 


1.441 


^3) 


0.715(2) 


1.542(25) 





0.21 


0.0109(1) 


0.0449 


.3) 


1.272 


'A) 


0.552(3) 


^ A rw f a\ 

1.405(34) 





0.22 


0.0134(1) 


0.0548 


^32) 


1.086 


[?,) 


0.400(2) 


1.300(37) 





0.23 


0.0161(2) 


0.0681 




0.871 


[h) 


0.253(2) 


1.074(65) 


3 


0.18 


0.2174(3) 


0.0252 


.13) 


1.554 


[7) 


0.808(5) 


1 ^ f A a\ 

1.631(44) 


3 


0.19 


0.2201(2) 


0.0361 


J) 


1.376 


J) 


0.621(5) 


1.473(48) 


3 


0.2 


0.2247(3) 


0.0454 


^ -1 A\ 

.14) 


1.179 


'5) 


0.437(3) 


1.342(61) 


3.5 


0.175 


0.2587(3) 


0.0268 


.12) 


1.606 


'5) 


0.848(5) 


"1 f\ ( c\ a \ 

1.679(24) 


3.5 


0.185 


0.2624(3) 


0.0333 


.13) 


1.408 


'7) 


r\ A r\ ^ A\ 

0.649(4) 


1.532(31) 


3.5 


0.195 


0.267(4) 


0.0436 


.13) 


1.211 


J) 


0.461(4) 


1.392(43) 


4 


0.17 


0.3034(1) 


0.0249 


.2) 


1.623 


[3) 


0.874(2) 


1.688(5) 


4 


0.18 


0.3079(1) 


0.0318 


.2) 


1.426 


[3) 


0.659(2) 


1.523(6) 


4 


0.19 


0.3141(1) 


0.0408 


.2) 


1.207 


^ A \ 

.4) 


0.458(2) 


1.367(14) 


4 


0.2226 


r\ A r\r\c\ ( '—!\ 

0.4002(7) 


0.1328 


^24) 


0.831 


.35) 


-.074(11) 




A -1 

4.1 


0.2211 


r\ A ^ r\r\ ( A \ 

0.4300(4) 


0.1345 


^20) 


0.997 


.41) 


-.100(9) 


1 r\''y ( A a\ 

1.07(44) 


4.2 


0.2195 


r\ A A A ^ / /^\ 

0.4445(6) 


0.1535 


'29) 


1.254 


'32) 


-.081(9) 


2.12(36) 


4.3 


0.165 


0.3319(2) 


0.0220 


J) 


1.663 


J) 


r\ r\c\r\ ( r'\ 

0.920(5) 


1.715(21) 


4.3 


0.175 


0.3367(2) 


0.0293 


;8) 


1.463 


[7) 


0.696(5) 


1.546(22) 


4.3 


0.185 


0.3440(2) 


0.0404 


;8) 


1.242 


:6) 


0.485(3) 


1.379(24) 


4.3 


0.205 


0.3732(3) 


0.0736 


;8) 


0.647 


:8) 


0.094(3) 




4.3 


0.21 


0.4286(4) 


0.1369 


;i2) 


0.755 


'74) 


-.026(9) 




4.3 


0.218 


0.4661(4) 


0.1790 


;i8) 


1.413 


;i3) 


-.083(9) 


1.784(55) 


4.5 


0.16 


0.3524(2) 


0.0209 




1.732 


;6) 


0.997(4) 


1.782(8) 


4.5 


0.17 


0.3580(2) 


0.0282 


;6) 


1.520 


;6) 


0.760(4) 


1.595(1) 


4.5 


0.18 


0.3656(2) 


0.0384 




1.298 


;5) 


0.534(4) 


1.423(8) 


4.5 


0.195 


0.3856(3) 


0.0590 


;9) 


0.882 


:ii) 


0.201(4) 


1.145(38) 


4.5 


0.2 


0.4007(6) 


0.0807 


19) 


0.696 


;24) 


0.090(5) 





Table 13: Results of the plaquette, the Polyakov loop, the pion screening 
mass, twice the quark mass, and the rho meson screening mass for = 2 
obtained on an 8^ x 10 x 4 lattice. Data marked with f are taken from 
Ref. iT3i obtained on an 8^ x 20 x 4 lattice. 
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p 


K 


plaquette 


Polyak 


ov 






2m^a 


rUpa 


4.5 


O.zOz 


0.4591 


^3) 


0.1643 


.9) 


1 IOC 

1.135 


^36) 


-.072(8) 




4.0 


U.zUo 


0.4752 


^3) 


U.ioU9 




1 yl 01 

i.4zi 


.18) 


1 oo otN 
-.iz8(z/) 


1 I'OQ/ AO\ 

1. /3o(4o) 


4.0 


n o 1 /I o 
U.zi4o 


0.4949 




U.zio 1 


'1 «^ 
^i6 ) 


1 o 
i.OOZ 


'1 n\ 

.10) 


no /I /'oN 
-.U34(8 ) 


i. /4D(iy ) 


r 



n 1 /I 


U.4U95 


^zj 


U.Uiz8 


J) 


o n /I 
Z.U4o 


.5) 


1 o 7n f K\ 

i.3 /y(o j 


o n7o//;^N 
z.U i Z(0 ) 


r 



U.iO 


U.4148 




U.Uzi / 


J) 


1 oni 
i.8Ui 


.12) 


1 nno /'I N 

i.uyoi^ij 


1 ooo 1 nN 
i.8z8(iy j 


r 



U.iD 


0.4zl5 


^zj 


U.UoUi 


J) 


1 c;i 

i.ooi 


.8) 


n one; /'rN 
U.oUol^O ) 


1 /^n /I 1 N 
i.0U4(i j 


r 



n 1 T 


0.4351 


^zj 


U.U4Uy 




1 oTn 
i.z /y 


.6) 


n KOO^ciN 
U.ozz(^o ) 


1 QO/I ^ 1 ON 

i.oy4(iz ) 


r 



U.io 


0.5174 




U.zzoU 




1 /I on 
i.4oU 




-.U80(^y j 


1 CO^ 1 1 1 N 

i.bo / (ii) 


r 



u.iy 


U.5378 


^Ij 


U.Z08U 


J) 


i.686 


'1 1 \ 
.11) 


-.Uy6(8 ) 


1 O^^ 1 N 

i.86i(8 j 


r 



U.iyoz 


U.5473 


^1) 


n o'Ton 

U.z /8y 




1 Vl 7 


.5) 


n non f 
U.Uzy(^0 j 


i.8o / ( / j 


O.ZO 


U.i 


0.44zD 


^zj 


n nni q 
U.UUio 




o nQ/i 


.5) 


z.oU8(^0 ) 


Z.yo / (0 ) 


O.ZO 


nil 
U.ii 


r\ A A A a 
0.4446 


^zj 


n nn/i i 
U.UU4i 




O RQ1 
Z.DOi 


.10) 


o 1 ni /"oN 
z.iyi^Sj 


O <^Q7/'l 1 N 

z.Do/(ii) 


D.ZO 


n 1 o 
U. iz 


r\ A A '70 

U.447z 


^z) 


U.UU80 


J) 


/I OQ 
Z.4ZO 


.') 


1 Qfi7/'7N 
i.oD / / j 


O /I QQ /'fiN 


O.ZO 


U.io 


A /I TAO 

U.45Uz 


^zj 


n n 1 /I n 
U.Ui4U 


J) 


O 1 O /I 

z.i84 


.8) 


i.o63(^6 ) 


o onn/nN 
z.zUU(y ) 


O.ZO 


U. 14 


U.4556 


^ZJ 


n noi Q 
U.Uzio 


J) 


1 n /1 1 

i.y4i 


.4) 


1 ORQ /'/IN 

i.z0o(^4j 


1 n7n/' c^N 
i.y ( U(0 ) 


O.ZO 


U.io 


0.4635 


^3j 


n none; 
U.UoUo 


.9) 


i.DO/ 


'1 ON 

.12) 


n non/'nN 

u.yoy(^y ) 


1 7nn/' 1 /I N 
i. /Uy(i4j 


O.ZO 


U.ioo 


0.4746 


^3) 


n n/inn 

u.U4yy 


'1 o^ 


1 /I nc: 

i.4yo 


.6) 


U./OD(D) 


1 c:^;q^7N 
i.ODo( / ) 


O.ZO 


U.io 


A A O A 

0.4846 


^3) 


n n^^To 
U.UO / 8 


'1 1 ^ 


1 QO /I 

i.oz4 


'nN 
.9) 


U.O /U(^ / ) 


1 Qn7/'i oN 


O.ZO 


U.ioo 


A r O A 'V 

0.5307 


^z) 


n oo /1 1 
U.zz4i 


'1 n^ 
.10) 


i.3oi 


.9) 


n 1 7o /nN 
U.i (3(9) 


1 /I^^O/ION 

i.4o8(iz j 


O.ZO 


U.i / 


0.5513 




n Oi^nc^ 

u.zoyo 


.8) 


1 C^QI 

i.ooi 


'1 ON 

.12) 


1 OR /"I ON 

-.ic)D(^i8 ) 


1 ^^n^ z' 1 nN 
i.Dyo(iU j 


O.ZO 


n 1 o 
U.io 


U.oooy 


^1 J 


n ooi^ 1 
U.z56i 




i.6y6 


'RN 

.6) 


-.i6U(6) 


1 o o /^N 

i.8o3( / ) 


0.0 1 


U. iO 


0.5530 


[2) 




.') 




.0) 


n c:i on 


i.OZo(^o J 


5.5t 


0.16 


0.5662 


[2] 


0.2815 


;5) 


1.415 


7) 


0.103(5) 


1.490(7) 


5.5t 


0.1615 


0.5677 


[2) 


0.2863 


;8) 


1.441 


;5) 


0.048(6) 


1.513(9) 


5.5t 


0.163 


0.5699 


;i) 


0.2905 


;8) 


1.438 


;8) 


-.016(4) 


1.506(8) 


6 


0.15 


0.6122 


;2) 


0.3456 


;io) 


1.469 


;7) 


0.233(5) 


1.510(8) 


6 


0.1524 


0.6131 


;3) 


0.3478 


:i6) 


1.467 


7) 


0.142(8) 


1.514(7) 


6 


0.155 


0.6157 


[2) 


0.3555 


;9) 


1.480 


;5) 


0.042(7) 


1.529(9) 


6 


0.16 


0.6188 


[2) 


0.3607 


;9) 


1.534 


;6) 


-.120(6) 


1.594(8) 


10 


0.13 


0.7853 


;i) 


0.6126 


;ii) 


1.496 


;6) 


0.447(2) 


1.491(6) 


10 


0.14 


0.7865 


;i) 


0.6157 


;8) 


1.439 


;8) 


-.010(4) 


1.437(9) 


10 


0.15 


0.7873 


;i) 


0.6230 


:8) 


1.591 


:2) 


-.427(6) 


1.598(3) 



Table 13: Continued. 
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K 


plaquette 


Polyakov 




2m^a 


TUpa 


4.2 


0.2195 


0.4410(3) 


0.0067(12) 


0.897(50) 


-.082(11) 






4.3 


0.2183 


0.4593(2) 


0.0054(6) 


1.070(62) 


-.111(21) 






4.4 


0.2163 


0.4742(1) 


0.0098(7) 


1.241(28) 


-.105(8) 


1 


677(52) 


4.5 


0.2143 


0.4889(2) 


0.0071(10) 


1.371(11) 


-.058(9) 


1 


604(31) 


5 


0.1982 


0.5455(1) 


0.0831(10) 


1.638(7) 


0.044(9) 


1 


749(8) 


5.02 


0.16 


0.4256(1) 


0.0023(4) 


1.542(5) 


0.800(4) 


1 


605(6) 


5.02 


0.17 


0.4384(1) 


0.0041(3) 


1.242(6) 


0.508(4) 


1 


343(11) 


5.02 


0.18c 


0.4696(1) 


0.0102(5) 


0.710(10) 


0.149(5) 





986(29) 


5.02 


0.18d 


0.5180(2) 


0.0399(9) 


0.923(47) 


-.164(18) 


1 


40(16) 



Table 14: The same as Table ^ for Np = 2 on a 12'^ x 6 lattice. 



K 



plaquette Polyakov 



2mga 



nipa 



4.5 0.16 0.3522(1) 0.0010(7) 1.731(7) 0.999(5) 1.779(20) 

4.5 0.17 0.3574(1) 0.0004(6) 1.513(5) 0.759(4) 1.588(19) 

4.5 0.18 0.3649(1) 0.0023(6) 1.281(5) 0.529(3) 1.398(18) 

5.5t 0.15 0.5377(3) 0.0073(1) 1.115(16) 0.542(6) 1.167(19) 

5.5t 0.155 0.5481(3) 0.0081(2) 0.807(35) 0.308(14) 0.874(39) 

6 0.1524 0.6131(3) 0.3478(16) 0.837(19) -.003(4) 0.881(22) 



Table 15: The same as Table |T3| for Np = 2 on an 8^ x 10 lattice. Data 
marked with f are taken from Ref. [|l^ obtained on an 8^ x 20 lattice. 
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0.2226 


0.4944 
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0.2637 
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0.009(4) 


1 


689(6) 


4.5 


0.16 


0.3598 
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0.0334 


7) 
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0.3691 


[2) 


0.0459 


;6) 
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4.5 
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[2) 
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;7) 
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0.478(4) 


1 


385(14) 


4.5 
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[2) 


0.0812 
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1 
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1 
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Table 16: The same as Table |T3| for Np = 3 on an 8^ x 10 x 4 lattice. 
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Table 16: Continued. 
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Table 17: The same as Table O for Np = 3 on a 12'^ x 4 lattice 
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Table 18: The same as Table O for Np = 6 on an 8'^ x 10 x 4 lattice. 
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Table 19: Results of the plaquette, the Polyakov loop, the pion screening 



mass, the rho meson screening mass, and the phi meson screening mass for 
Np = 2 + 1 obtained on an 8^ x 10 x 4 lattice. 
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Table 20: The same as Table 19 for Np = 2 + 1 obtained on a 12 x 4 lattice. 
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Table 21: The same as Table IT^ for Np = 2 + 1 obtained on an 8 x 10 lattice. 
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Table 22: Hopping parameters for Np = 2 + 1 simulations performed on 

8^ X 10 X 4 and 12^ x 4 lattices. Kud for u and d quarks is set to be equal to 
Kc and Kg for s quark is chosen so that ~ 150 MeV and 400 MeV in the 
left and right columns, respectively. 
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Figure 1: Squared pion screening mass m^a^ and twice the quark mass IrHqa 
for = 2 at /3 = on an 8^ X 10 X 4 lattice. Errors are smaller than the 
size of symbols. Solid curves are the results of a strong coupling calculation, 
Eq.(@). 
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Figure 2: Schematic graph for the chiral hmit hne Kc defined by = or 
m^r = at T = in the coupUng parameter space {P, K). Also plotted are 
the curves where ruq = and m.,^ = at finite Nt- See text for discussions. 
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Figure 3: Physical quantities for Np = 2 at /3 = 5.0 on an 8^ x 10 x 4 lattice 
(open squares): (a) m^a^ and 2mqa, (b) the plaquette and the Polyakov loop. 
Plotted together are the data by the MILC collaboration on an 8^ x 20 x A^^ 
lattice with Nt = 4 (filled squares), 6 (triangles), and 8 (diamonds) [|1^]. 
The finite temperature transition Kt obtained by the MILC data locates at 
K = 0.177 — 0.178 {1/K = 5.62 — 5.65) for A^^ = 4. 
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Figure 4: The same as Fig. ^ at /3 = 
Nt = 4: (squares), 6 (triangles), and 8 
transition Kf locates at K = 0.200 — 

51 



4.5 on 8^ X 10 X A^^ lattices, where 
(diamonds). The finite temperature 
.202 {1/K = 4.95 — 5.0) for Nt = 4. 



0.24 I 




Figure 5: Phase diagram for Np — 2. Filled symbols are for Kc determined 
by ttItt — and rUq — 0. Open symbols are for Kt{Nt — 4) and other symbols 
such as crosses except filled ones are for Kt{Nt = 6). Circles are our data. 
Lines are to guide the eyes. 
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Figure 6: Molecular-dynamics time history of A^inv for Np = 2 on the line 
obtained on an 8^ X 10 X 4 lattice. 
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Figure 7: The pion screening mass squared m^a for Np = 2 on the Kc hne 
obtained on 8^ x 10 x 4 and 12^ x 6 lattices. 
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Figure 8: Phase diagram for Np — 2 and 3. Filled symbols are for Kc{m'^) 
and Kc{mq). Open symbols are for Kt{Nt = 4) for A^^ = 2 and open circles 
with cross for Np = 3. Circles arc our data. On the Kt line for Np = 3, clear 
two-state signals are observed at P < 4.7 both on 8^ x 10 x 4 and 12^ x 4 
lattices. Lines are to guide the eyes. 
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Figure 9: Time history of A^inv for Np = 3 on the Kc hne obtained on an 
8^ X 10 X 4 lattice. 
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Figure 10: m^a^ for Np = 3 on the K^. line obtained on an 8^ x 10 x 4 lattice. 
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Figure 11: Time history of A^inv for Np = 6 on the K^. line obtained on an 
8^ X 10 X 4 lattice. 
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Figure 12: m^a^ for Np = 6 on the Kc line obtained on an 8^ x 10 x 4 lattice. 



59 




Figure 13: Physical quantities for Np = 3 at /5 = 4.0 on 8^ x 10 x 4 (open 
squares) and 12'^ x 4 (filled diamonds) lattices: (a) m^a^ and 2mqa, (b) 
the plaquette and the Polyakov loop. The finite temperature transition Kt 
locates at K ~ 0.204 {1/K ~ 4.90). 
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Figure 15: The same as Fig. O at /3 = 4.7 obtained on a 12^ X 4 lattice. The 
finite temperature transition locates at i^T ~ 0.1795 {l/K ~ 5.57). 
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Figure 16: The same as Fig. |1^ at /5 = 5.0: (a) m^a^ and 2mga, (b) the 
Polyakov loop. The finite temperature crossover Kt locates at ~ 0.166 — 
0.1665 {1/K ~ 6.01 — 6.02). 
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Figure 17: The same as Fig. O at /5 = 5.5 obtained on an 8^ x 10 x 4 
lattice. The finite temperature crossover Kt locates at K ~ 0.125 — 0.130 
(l/ir~ 7.7 — 8.0). 
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Figure 18: Time history of the plaquette for Np = 3 at (a) (3 = 4.7 and 
5.0 on a 12^ x 4 lattice. 
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Figure 19: The rho meson screening mass m^a and twice the quark mass 
2mqa in the confining phase as a function of — l/Kc- Open symbols are 
for Nf ^2, p ^ 3.0, 3.5, 4.0, 4.3, and 4.5 on an 8^ x 10 x 4 lattice. Filled 
symbols are for AT^ = 3, /3 = 4.0, 4.5 and 4.7 on 8^ x 10 x 4 and 12^ x 4 
lattices. The values of Kc{(3) for Np = 2 is used. Horizontal errors are from 
those for Kc with taking into account the difference due to definitions, either 
the vanishing point of or nig- 
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Figure 20: Time history of (a) A^inv and (b) the plaquette for ~ 150 MeV 
on an 8^ X 10 X 4 lattice. 
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Figure 21: Time history of the plaquette for ~ 400 MeV on a 12^ x 4 
lattice. 
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Figure 22: m^a^ versus /3 for ~ 0, 150 and 400 MeV with rUud — 0. Filled 
and open symbols are for 8^ x 10 x 4 and 12^ x 4 lattices, respectively. 
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Figure 23: Order of the finite temperature QCD transition in the {mud-,TT^s) 
plane. First order signals are observed at the points marked with filled circle, 
while no clear two state signals are found at the points with open circle. The 



second order transition line is suggested to deviate from the vertical axis 
as rriud oc (m* — m^)^/^ below m*. The values of quark mass in physical units 
are computed using determined from nip. a^^ ~ 0.8 GeV for (3 < 4.7 and 
~ 1.0(1.8) GeV for (3 = 5.0(5.5). See Sec. ^ for more detailed discussion on 
the values of the quark mass in physical units. The real world determined by 
the value of m^/rrip and nin/mp corresponds to the point marked with star. 
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